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Nonequilibrium Kondo transport through a quantum dot in a magnetic field 
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We analyze universal transport properties of a strongly interacting quantum dot in the Kondo 
regime when the quantum dot is placed in an external magnetic field. The quantum dot is described 
by the asymmetric Anderson model with the spin degeneracy removed by the magnetic field resulting 
in the Zeeman splitting. Using an analytical expression for the tunneling density of states found 
from a Keldysh effective field theory, we obtain in the whole energy range the universal differential 
conductance and analytically demonstrate its Fermi-liquid and logarithmic behavior at low- and 
high- energies, respectively, as a function of the magnetic field. We also show results on the zero 
temperature differential conductance as a function of the bias voltage at different magnetic fields as 
well as results on finite temperature effects out of equilibrium and at a finite magnetic field. The 
modern nonequilibrium experimental issues of the critical magnetic field, at which the zero bias 
maximum of the differential conductance starts to split into two maxima, as well as the distance 
between these maxima as a function of the magnetic field are also addressed. 

PACS numbers: 73.63.Kv, 72.10.Fk, 72.15.Qm 
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I. INTRODUCTION 

Among various quantum many-particle phenomena in 
condensed matter physics the Kondo effect^ stands out as 
a unique fundamental state exhibiting remarkable prop- 
erties both in bulk systems^ii^ and in a quantum dot (QD) 
setup^""— . 

In QDs^ the Kondo effect appears when a QD is 
switched into the Coulomb blockade regime by tuning 
the gate voltage. Normally the QD states with even and 
odd numbers of electrons (Coulomb valleys) result in a 
zero bias minimum of the differential conductance as a 
function of the bias voltage. However, when the tempera- 
ture T is lowered this situation changes qualitatively: the 
QD states with even numbers of electrons still result in a 
zero bias minimum of the differential conductance while 
for the QD states with odd numbers one observes an 
anomalous behavior where the differential conductance 
at zero bias shows a maximum often referred to as the 
Kondo resonance. This so-called zero bias anomaly is a 
manifestation of the unpaired spin inherent to the QD 
states with odd numbers of electrons. 

Another important property of the Kondo effect is that 
together with its appearance the zero bias anomaly ac- 
quires universality: the differential conductance of differ- 
ent QDs collapses into a single function of the temper- 
ature and voltage. It depends on specific QD parame- 
ters only through the energy scale fcTx, where Tk is the 
Kondo temperature. 

The Kondo effect has different universal character at 
low and high energies. At low energies, smaller than fcTx, 
one observes the Fermi-liquid behavior (strong coupling 
limit), a quadratic universal dependence of the differen- 
tial conductance on the temperature, voltage, or external 
magnetic field. In the opposite regime of high energies, 
larger than fcTx, the differential conductance is charac- 
terized by a logarithmic universal dependence (weak cou- 
pling limit). At intermediate energies it shows a universal 



crossover between these two limits. 

The zero bias anomaly and its universality have been 
observed in equilibrium- and nonequilibriunii^. The 
nonequilibrium Kondo universality in the deep crossover 
is characterized by a universal value of the differential 
conductance at the Kondo voltage, eV — kTa- Theory 
predictsiiii^ that this universal value is equal to 2/3 of 
the maximum value which the differential conductance 
takes in equilibrium at zero temperature. This provides 
an efficient way to nonequilibrium experimental measure- 
ments of the Kondo temperaturoH. 

Interplay between the Kondo state and other quan- 
tum collective phenomena such as ferromagnetismM, 
superconductivityi^^ or the Kondo state in an exter- 
nal magnetic fiel d^ '^^ currently represent a very active 
research field both in theory and experiment. 

The Kondo effect in a QD placed in an external mag- 
netic field is of particular interest. The magnetic field 
couples to the spins of the electrons in the QD and pro- 
duces the Zeeman splitting of the single-particle QD en- 
ergy level. It is known^"— that such a magnetic field 
destroys the zero bias anomaly, that is, it restores the 
normal behavior where the QD states with odd numbers 
of electrons at low voltages result in minima of the dif- 
ferential conductance, similar to the QD states with even 
occupancies. 

However, it turns out that the magnetic field does 
not destroy the universality. Theorji and experiments^ 
demonstrate that even when there appears a minimum in 
the differential conductance at 1/ = 0, it still remains uni- 
versal both in equilibrium (as a function oi T at V = 0) 
and nonequilibrium (as a function of T and V). 

Important nonequilibrium issues which have been ad- 
dressed in experimentsiiii^ on the Kondo effect in an 
external magnetic field are 1) Kondo universality; 2) the 
critical magnetic field at which the zero bias maximum of 
the differential conductance starts to split into two max- 
ima; 3) the distance between these maxima as a function 



of the magnetic field; 4) the high-field limit of this dis- 
tance. The issues above have obviously a highly nonequi- 
librium nature. The nonequilibrium theory of the Kondo 
effect in the whole energy range is a complicated prob- 
lem even without a magnetic field. The numerical renor- 
malization group (NRG) method^- which is known to be 
a numerically exact tool in equilibrium cannot be easily 
generalized to nonequilibrium. To cover the whole energy 
range in nonequilibrium one has to resort to other meth- 
ods. These are the non-crossing approximation (NCA) 
(intermediate and large energies) ^'^°i^^ , equations of mo- 
tion (mainly qualitative tool)^ '^^i^^ , mean-field theories 
(low energies )2ii2^, or 1/7V expansions (low energies)^. 
In the whole energy range the real-time renormalization 
group methodic is a powerful tool. It has been devel- 
oped for the s — d modeli. To our knowledge, however, 
a generalization to include the Zeeman splitting has not 
been done yet. Another theoretical tool is the Keldysh 
effective field theorj*i^i22i^ developed for the Anderson 
model^. In its first application it has been shown that 
this theory well captures the weak coupling limit^ and 
is capable to cope with finite, though large interaction 
strengths™. Recently it has been further developed to 
describe the behavior of the differential conductance in 
the whole range of temperatures and bias voltages^S. 

In the present work we address the behavior of the 
Kondo state in an external magnetic field both in equi- 
librium and nonequilibrium. To this end we extend the 
Keldysh effective field theory^^ to account for the Zee- 
man splitting of the single-particle QD energy level and 
derive an approximate analytical expression for the QD 
tunneling density of states (TDOS). The theory is then 
applied to calculate the universal differential conductance 
as a function of the temperature, bias voltage and mag- 
netic field in the whole energy range. We demonstrate 
that despite the approximations made in the theoretical 
derivation the theory has a number of advantages: 1) 
Fermi-liquid behavior when the Zeeman energy is much 
less than /cTk (advantage over NCAi); 2) logarithmic be- 
havior when the Zeeman energy is much larger than fcTk 
(advantage over mean- field theories^); 3) critical mag- 
netic field in a good agreement with previous theories^S 
and experiments^^; 4) distance between the maxima of 
the differential conductance as a function of the mag- 
netic field in a good agreement with experimentsii; 5) 
high-field limit of this distance in a good agreement with 
experimentsii; 6) universality with the correct scaling 
given by Tk- 

The paper is organized as follows. In Section |TT] we 
present our model of a strongly correlated QD in a mag- 
netic field. Here we also show the slave-bosonic trans- 
formation used in our work. The analytical expression 
for the QD TDOS resulting from an effective Keldysh 
field theory is given in Section Hill Here we analytically 
investigate the low-energy sector of the theory proving 
its Fermi-liquid nature and providing the Fermi-liquid 
coefficients for the temperature, voltage and magnetic 
field behaviors as well as their universal ratios. Addi- 



tionally, when the Zeeman energy is much larger than 
/cTk we analytically derive the logarithmic asymptotics 
of the differential conductance as a function of the mag- 
netic field at r = and V = 0. In Section lTVl integrating 
the analytical expression for the QD TDOS, we calculate 
the universal behavior of the differential conductance in 
the whole range of temperatures, voltages and magnetic 
fields. Here we show all the regimes of the Kondo state: 
the strong coupling limit, weak coupling limit and the 
crossover region. Finally, we conclude the paper in Sec- 
tion |V] where all the advantages and drawbacks of our 
theory are summarized and a systematic way to improve 
the theory is mentioned. 



II. THEORETICAL MODEL 

To describe the strongly interacting QD we employ the 
single- impurity Anderson model22. (SIAM). This model 
describes a system of two interacting electrons. The 
strength of the electron-electron interaction is given by 
the parameter U > 0. In the absence of the interaction 
the electrons can occupy a spin degenerate single-particle 
energy level Cd. We include the effect of an external 
magnetic field through the Zeeman splitting of this en- 
ergy level, Ed — ^ Co- with eg. = ed + o'Ae/2, a = ±1, 
Ae = gfiBH, where g is the ^-factor, /ib is the Bohr mag- 
neton and H is the magnetic field. The QD Hamiltonian 
is thus of the following form: 



-f^QD — 



E- 



Ufifhi, 



(1) 



where the number operator fiu is given in terms of the 
original Anderson fermionic operators as 

fia^dldcy, {da,dl,} ^ da,a' {da,da'} = 0. (2) 

In addition to the interacting electronic system de- 
scribed by the SIAM, we also consider a noninteracting 
electronic system which models the contacts, that is an 
external system used to probe the QD. We consider two 
contacts and the corresponding Hamiltonian is 

k,a,x (3) 

X^kaxT ^f^'ff' x' f — ^kax.k' a' x' 7 X^kax^ ^k' a'x' j — U, 

where k is the set of orbital quantum numbers in the 
contacts, a is the contact spin degree of freedom and 
x = L, R denotes the left and right contacts. Here we 
also assume that the electronic states in the left and right 
contacts are characterized by the same complete set of 
quantum numbers, {k, a}. The left and right contacts are 
in equilibrium states specified by the chemical potentials 
/iL and /j,R,, respectively. 

The electrons can tunnel between the QD and con- 
tacts. These tunneling events are described in terms of 



the tunneling Hamiltonian, 

Ht^Y. (Tkaclja + nAcka.), (4) 

where Tk^ are the tunnehng matrix elements and we as- 
sume that these matrix elements are diagonal in the spin 
space and that they do not depend on x, that is, they are 
the same for the left and right contacts (QD symmetri- 
cally coupled to the contacts). 

In equilibrium the chemical potentials of the QD and 
contacts are equal to each other. We denote this equilib- 
rium value as /xg and define the applied bias voltage V 
as ^L,R = fJ-oT eV/2. 

As mentioned in the introduction, the Kondo effect, 
or the zero bias anomaly, develops at low temperatures 
out of the Coulomb blockade regime when the QD is 
blockaded with an odd number of electrons. Within the 
SIAM this means that the Kondo resonance takes place 
when the QD has one electron. This is achieved when 
the electron-electron interaction U exceeds the energy F 
(see below) resulting from the QD-contacts coupling. 

To describe the relevant physics and at the same time 
to simplify the formalism, we consider the case of of the 
strongly interacting QD where [/ — > oo. In this case 
the double occupancy is forbidden and the QD may have 
zero or one electron. At the same time the difference 
fiQ — Ed is assumed to be finite but much larger than F so 
that the QD has one electron. This model admits the so- 
called slave-bosonic representation^!.2ii^. It represents a 
transformation where the original QD fermionic opera- 
tors da, d\ are mapped onto new fermionic operators /a-, 
fl and bosonic (often called slave-bosonic) operators 6, 
b^ by means of the following relations: 



Ft 






{TkAaJ^b'^ +njlbcu..)- 



(8) 



xJcTi J a' > "(y.(J' 1 xJcTi let' I 

[6,6t] = l, [6,6]=0. 



(5) 



Additionally, the operators /^ and /t commute with the 
operators 6, h^ . 

Physically the transformation ([S]) means that instead 
of using the states of the electrons in the QD one uses 
the states of the QD itself: the empty state (&, W) and 
the state with one electron (/o-, fl). 

Since the empty state and the state with one electron 
represent all the states of the QD, one arrives to the 
constraint, 



flf, + h%^i, 



(6) 



which physically means that the total number of the new 
fermions and slave-bosons is restricted to be equal to one. 
After the slave-bosonic transformation the QD and 
tunneling Hamiltonians become: 



H, 



QD 



/ , ^aJaJin 



(7) 



The total Hamiltonian of this, in general nonequilib- 
rium, problem is the sum of Hqd, Hq and Ht given 
by Eqs. ([7]), @ and ([5]), respectively. This Hamiltonian 
together with the constraint in Eq. © represent the the- 
oretical model able to describe the essential behavior of 
the zero bias anomaly arising due to the Kondo effect in 
the presence of an external magnetic field. 



III. KELDYSH FIELD INTEGRAL SOLUTION 
AND ITS ASYMPTOTICS 

As it is well known, the Kondo effect is an essentially 
nonperturbative phenomenoni. Therefore, to cover all 
the regimes of the Kondo state one must resort to non- 
perturbative methods able to deal with both electron- 
electron interactions and nonequilibrium. 

In particular, the real-time formalism developed by 
Keldysb^^ represents a general and powerful method to 
treat interacting many-particle systems both in equilib- 
rium and nonequilibrium. Originally this method was 
developed in the diagrammatic framework. Its field in- 
tegral form has been developed in Ref. [SJ to study the 
interplay between disorder and electron-electron inter- 
actions in metals. The advantage of the field integral 
form is that it reduces the problem to an analysis of a 
functional, the so-called Keldysh effective action. This 
analysis has a more systematic nonperturbative charac- 
ter than the diagrammatic selection on the basis of the 
topological structure reflecting the physical content of a 
given diagram. 

In addition to being systematic, this nonperturbative 
approach turns out to be very general and applicable 
also to mesoscopic and QD systems. It has been applied, 
e.g., to describe the Coulomb blockade in QD a35,36 g^g 
well as the Kondo effect in the strong coupling limit^^, 
weak coupling limi t^^i^^ and in the whole energy rangei^. 

In particular, the slave-bosonic formulation presented 
in the previous section has been used in the Keldysh 
effective action framework in Refs. [l3l27ll28l The ad- 
vantage of these slave-bosonic theories over slave-boson 
mean-field theorie o"'^'^^ and 1/N expansions'^ is that they 
take into account the constraint in Eq. ^ exactly while 
slave-boson mean- field theories and 1/A^ expansions ac- 
count for this constraint only approximately. Because 
of this approximation one can only access the low energy 
(smaller than fcTx) physics of the Kondo state, that is the 
Fermi-liquid regime, or strong coupling limit, character- 
ized by the quadratic dependence of the QD differential 
conductance. The main drawback of this approximation 
is that it does not lead to logarithmic termsi in the differ- 
ential conductance and thus one cannot access intermedi- 
ate (of order of fcTx) and high energy (larger than kT^) 
physics, that is, the crossover region and weak coupling 
limit of the Kondo state, respectively. In contrast, the 
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the QD is allowed to have arbitrary numbers of the new 
fermions and slave-bosons, that is one deals with the full 
Hilbert space. One projects onto the physical QD Fock 
space by taking the limit /i — >■ oo. 

In Eq. ^ So[x''^''iity,x''^'Ht)] is the standard 
free bosonic action^^ on the Keldysh contour and 
ST[x''^'%t);x''^''^{t)] = -ihtT\n[-iG^°y-^-iT] is the tun- 
neling action of the problem. Here the matrix T is off- 
diagonal in the QD-contacts space, 



T = 





MT{kat\a't') 



Mr^{at\k'a't') 



(11) 



FIG. 1: (Color online) To formulate the field integral in Eq. 
(|10|l the constraint in Eq. (|6} is temporarily removed by 
introducing a positive real parameter /i. When /i is finite, 
< /^ < oo, the QD may have any number of the new fermions 
and slave-bosons. One projects onto the physical subspace, 
the QD Fock space, by taking the limit /i — > oo after the field 
integration has been performed. 



slave-bosonic Keldysh effective action theories in Refs. 
[l2,27.28 account for the constraint in Eq. © exactly 
and thus contain logarithmic terms. 

As shown below, the formalism of Ref. 12: is straight- 
forwardly generalized to the case of a QD in an external 
magnetic field. 

The theory of Ref. \l^ introduces a complex function 
Ea which is found from the condition that the QD TDOS, 



^^(e)^-^MG+.(e)], 
nil 



(9) 



where G^^, (e) is the QD retarded Green's function, in 
equilibrium and at zero temperature has a peak at e = /ig 
and the amplitude of this peak is equal to the unitary 
limit. 

In the case of a QD in an external magnetic field, one 
obtains the Keldysh effective action SeS using a deriva- 
tion similar to the one in Ref. |l2| and expresses any 
physical observable O = J^{d\, da) as 

(0)(t) = _L lim e'3^/-D[x,x]e^^"«l^""(*)'^""(*)lx 

J\Iq /i-J-OO J 

x^r'^^Wix^'-'^w], (10) 

+ 5T[x=^'^(i);x^'''^W], 

where x'^^''^ are the classical and quantum^S- eigenstates of 
the bosonic annihilation operator b, /3 = 1/kT is the in- 
verse temperature and A/q is a normalization constant^^. 
The limit /i — > oo in Eq. pH)) takes into account the 
constraint in Eq. ([5]). The point is that the field integral 
in Eq. (jlOp is most simply formulated for a system with 
arbitrary numbers of the new fermions and slave-bosons, 
that is temporarily removing the constraint in Eq. ([5]). 
One achieves such a formulation introducing a real pos- 
itive parameter fi. As shown in Fig. [U for < /i < oo 



MT{k(Tt\a't') 



S{t - t')5aa,Tko 



r\t) - 7,72 

x'^it) 



X' 



V2h 

xHt) 

^'(t)-7.V2. 



(12) 



The Green's function matrix G'"'' is block-diagonal in the 
QD-contacts space. Its QD block GJj has the standard 
2x2 fermionic Keldysh structure: 



Gf\at\a't') 



G+(i-i') GK(i-i') 
G-it-t') 



(13) 



In the frequency domain the components of the above 
matrix are 



G-{lu) = [G+i^W, 

G^H = J[G+(a.)-G-(c.)]^tanh 



(14) 



huj — /ij 



2kT 



Here Ea,a = ao-F/2, aa = "f^Sa, F = 27ri'c|rp {vq is the 
contacts density of states, r is the value of the tunneling 
matrix element Tka assumed to be independent of /c, a). 
In the case of the QD TDOS for spin a the expression 
for the integrand in Eq. (ITU)) is 



J-r''''(f);x"''^(i)] = [X-(t)x+(0) 
x[G(o)-i+r]-i(at|aO), 



X+(t)x-(0)]x 



(15) 



where x± ? X± are the slave-bosonic fields on the forward 

and backward branches^^ of the Keldysh contour. 

Anticipating that the Kondo physics arises from the 

spin-flip processes, to obtain the QD TDOS for fixed 

spin a within the lowest order expansion (see below), 

we first choose in the a- and — cr-blocks of the ma- 

trix g(°) 
I 



as 



EL. = E. 



Tj the imaginary parts of Ea a and Ea^- 

^}a — crAe and E\_^ — E\, while the real 
parts are E^^^ = ^a.-a = E^- Here Ea is found, as in 
Ref. IT2I, at zero magnetic field, yielding E^/kTx ~ 1.42, 
Ei/kTK w 1.05 with kTK = 2Wexp[-7r{no - ed)/r] (W 
is the Lorentzian width of the contacts density of states) . 

After this choice the a- and —cr-blocks of GJj are both 



expressed through the spin —a. Finally, as in Ref. Il2l . 
we expand S'eff and T up to the second order in the slave- 
bosonic fields in x'^K*) - ^<y^, X^K*) -laV^ and x'^it), 

xHt)- 

The resulting QD TDOS obtained after performing the 
Gaussian field integral in Eq. ((TU)) is then given by an 
expression similar to the one in Ref. [l^ except for the 
spin dependence arising due to the Zeeman splitting, 



Vaie) 



(16) 



2n [e^ + aAe/2 e + T^+^ie^ + [rS+ (e)]^ ' 






2tt^ 



"1 
2 



^ 



W 



2 2TTkT 



All 



^ 



W 
27rfcT 



i(e-o-Ae) 



27rfcr 27rfcT 



1 



27rfcT 



(17) 



exp( 



kT 



+ 1 



where ijj is the digamma function. 

Here we would like to mention the main drawback of 
the approximations made in deriving Eqs. ([TB)) and ([TT]) . 
Due to the second order expansion of Eq. ([15]), the theory 
does not take into account all inelastic cotunneling pro- 
cesses which are of the fourth order in the tunneling ma- 
trix elements. These processes are included only partly, 
through the effective action. Therefore, the present the- 
ory underestimates ^'t(e) at energies e — /zq > gfJ^BH 
{H > 0) and ^'^.(e) at energies e — fiQ < —gfiBH. Further- 
more, it is known^^'^* that for the s — d model at very 
smaU magnetic fields igfJ,B\H\ < fcTx) the QD TDOS 
peak is located at e = {2/3)g^BH and at e = .g^B-ff at 
larger fields. Our theory for the highly asymmetric SIAM 
predicts, as one can see from Eqs. P^ and ([TT]), that the 
peak is always located close to e = g^^H . This perhaps 
can be attributed partly to the quality of our approxi- 
mation and partly to the physical difference between the 
s — d model and the highly asymmetric SIAM. It is tech- 
nically complicated to include higher order terms in Eq. 
(IT51) . and the corresponding theory will be the focus of 
our future research. However, below we demonstrate that 
already this simple theory has a number of advantages 
and is in a good agreement with some other theoretical 
predictions as well as with experiments. 

With the QD TDOS ^ we can calculate the differen- 
tial conductance using the expression^*' for the current 
through the QD, 



/oo 
de[nR(e)-nL(e)] 
-OO 



W^ 



4e^ 



w- 



:'^a[e), 



'^l.rXc) 



1 



(18) 



exp[/3(e - /io ± eV/2)] + 1 ' 



The differential conductance ad is obtained from Eq. 
T5)) via the derivative of the current through the QD 
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FIG. 2: (Color online) The QD universal differential conduc- 
tance as a function of the magnetic field in the Kondo regime 
at T = and V = 0. We use ^o - td = TT, VK = lOOF which 
gives /cTk/F « 5.63 • 10"*. The Kondo temperature Tk is 
defined in the text. The black solid curve shows the differen- 
tial conductance obtained from Eqs. (|16p . (|17p and (|18p . The 
dashed lines show the Fermi-liquid (red/upper) and logarith- 
mic (blue/lower) asymptotics given by Eqs. (I19p and (|22ll . 
respectively. 



with respect to the bias voltage, ad = dl/dV. The dif- 
ferential conductance is in general a function of the tem- 
perature, voltage and magnetic field, ad — ad(T,V, H). 
It can be obtained substituting the QD TDOS, Eq. (fT6|) . 
into the expression for the current through the QD, Eq. 
(dH). This is a difficuff task because the QD TDOS dill) is 
already a complicated function which, in addition, must 
be integrated. Therefore, in general one has to perform 
a numerical integration to obtain the differential conduc- 
tance. 

However, at low energies it turns out to be possible to 
obtain the differential conductance analytically. From 
Eqs. dUl), (HH) and (HH) it follows that at T < Tr, 
e\V\ < /cTk and gfj,B\H\ < fcTx the differential conduc- 
tance takes the form: 



adiT,V,H) = 



2e' 



i-oH^)' 



'Cv 



eV 

mi 



^cJ'J^\ 



2n 



V kTK ) 



(19) 



which shows that at low energies the QD TDOS (jT6|) 
leads to two important consequences. First, the differ- 
ential conductance is a universal function of the tem- 
perature, bias voltage and magnetic field with the cor- 
rect scaling given by the Kondo temperature Tk of the 
corresponding problem without magnetic field. Second, 
the differential conductance does not contain terms lin- 
ear in T, V and ff, which means that the QD TDOS 
P^ correctly reflects the low energy physics leading to 
the Fermi-liquid behavior. 



The Fermi-liquid coefficients are given as 



4 21n(2|fQ 



1 1 41n(2|£„|) + l 



"' 3 |£„P ' ^^ .2 


l^aP 


4 ln(2|£,|) + l 




with the foUowing universal ratios 




cy _ 3 41n(2|f„|) + l ^ 
CT 27r2 41n(2|£„|) + 2 '" 


:^ 0.13043, 


CH 3 21n(2|£„|) + 2 ^ 
CT 27r2 21n(2|£:a|) + l ^ 


a 0.19509, 



(20) 



(21) 



where £"„ = Ea/kTK- For comparison we also give 
the universal ratios for the symmetric Anderson model, 
where one first puts fj-a — e^ = U/2 and only after that 
takes the limit U ^ oo: cy/cr = 3/(271^) w 0.15198, 
ch/ct = l/7r2 « 0.10132. 

In equilibrium and at zero temperature one obtains 
from Eqs. ([TE| . p7)) and (IT51) the differential conductance 
at gfiB\H\ > fcTx: 



ad(0,0,i/) 



/i 4 in2 



' ffMBl-gh 



(22) 



As one can see, at large magnetic fields the differential 
conductance at T = and y = is also a universal 
function of the magnetic field with the correct scaling. 
Moreover, it correctly reflects the high energy physics 
leading to the logarithmic behavior -'^. 

Already these low- and high-energy asymptotics 
demonstrate advantages of the present theory over other 
slave-bosonic theories like mean-field theories or NCA 
which do not have the logarithmic or Fermi-liquid be- 
havior, respectively^. 

In the next section we obtain the differential conduc- 
tance from Eqs. p6)l . (fTTl) and (fTSi) in the whole energy 
range both in equilibrium and nonequilibrium. However, 
to obtain it in this case one has, in general, to perform 
the integration in Eq. ([T5)) numerically. 



IV. UNIVERSAL RESULTS IN THE WHOLE 
ENERGY RANGE 
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FIG. 3: (Color online) The universal differential conductance 
as a function of the bias voltage at zero temperature and 
different magnetic fields. The parameters are the same as 
for Fig. [2l When the magnetic field increases the zero bias 
anomaly first decreases and then starts to split into two peaks. 
This happens at the critical magnetic field gfisHc ~ 0.9 fcTx. 
The inset shows the distance between these two peaks as a 
function of the magnetic field. The circles in the inset show 
the experimental data from Ref. [ITI . 




FIG. 4: (Color online) The universal differential conductance 
as a function of the bias voltage at finite temperatures and 
magnetic field gfisH = IO/cTr. The parameters are the same 
as for Fig. [21 The increase in temperature leads to two effects: 
1) the height of the peaks reduces; 2) the width of the peaks 
increases (broadening). 



It has been demonstrated in Ref. 12 that the differ- 
ential conductance resulting from the Keldysh effective 
action theory is a universal function of the temperature 
and bias voltage in the whole energy range. It has also 
been shown that it has the correct scaling given by the 
Kondo temperature Tk- Here we provide the differential 
conductance resulting from the Keldysh effective action 
theory taking into account an external magnetic field, 
Eqs. (HH) and (HH). Using Eq. dUl), we obtain the linear 
conductance as a function of the magnetic field at T = 0, 
ad{0,0, H). It is shown in Fig. [21 As one expects, in 



the Kondo regime ad{0,0, H) turns out to be a univer- 
sal function of g^BH/kTK (or g/iB^/fcTk, where Tk is 
defined as fTd(fK, 0,0) = (7^(0, 0, 0)/2, Tk ~ 1.47 Tk). 

Let us now address the nonequilibrium Kondo physics 
in the presence of an external magnetic field. In Fig. [3l 
we show the zero temperature differential conductance 
as a function of the bias voltage at different magnetic 
fields. Again all the curves are universal with the scal- 
ing given by the Kondo temperature, which means that 
Eqs. HH) and ([T7]) provide the differential conductance 



with the correct scahng also in nonequihbrium and finite 
magnetic fields. From Fig. [3] one observes that when 
the magnetic field is increased the zero bias Kondo peak 
first decreases and then splits into two peaks. Our the- 
ory predicts that the zero bias peak starts to split at the 
critical magnetic field q^bHc ~ 0.9 /cTk which is in a 
good agreement with another theoretical estimateSS, pre- 
dicting gfiBHc w 1.0 A: Tk. The inset of Fig. |3] shows the 
distance A between the two peaks as a function of the 
magnetic field. At low magnetic fields A has a rapid in- 
crease and then saturates at A/2g/iB-ff ~ 1.017. This 
high field limit of our theory is in a good agreement with 
the theoretical prediction made in Ref. |33. As one can 
see from the inset of Fig. |31 our theory is also in a good 
agreement with experiments in Ref. [ij. In contrast, the 
theories in Refs. ISOlSTI both predict a saturation at much 
larger fields than observed in Ref. [13. 

Finally, in Fig. 2] we show the differential conductance 
as a function of the bias voltage at finite temperatures 
and at the magnetic field value g/iB-ff — 10 fcTx. Also in 
this case the differential conductance is a universal func- 
tion with correct scaling given by the Kondo tempera- 
ture. As one expects, when the temperature increases 
the two peaks in the differential conductance are lowered 
and broadened. 



V. CONCLUSION 

In conclusion, we have applied the Keldysh effective 
action theory to study the Kondo effect in QDs in an ex- 
ternal magnetic field both in equilibrium and nonequihb- 
rium. To this end we have generalized the Keldysh effec- 
tive action theory from Ref. |Ij to take into account the 



effect of the magnetic field through the Zeeman splitting. 
An approximate analytical expression for the QD TDOS 
has been obtained. We have used this QD TDOS to cal- 
culate the differential conductance in the whole energy 
range and demonstrated that the theory has a number 
of advantages: 1) Fermi-liquid behavior when the Zee- 
man energy is much less than /cTk (advantage over the 
non-crossing approximation (NCA)-); 2) logarithmic be- 
havior when the Zeeman energy is much large than /cTk 
(advantage over mean- field theoriesi); 3) critical mag- 
netic field in a good agreement with previous theories^ 
and experimentsii; 4) distance between the maxima of 
the differential conductance as a function of the mag- 
netic field in a good agreement with experimentsii; 5) 
high-field limit of this distance in a good agreement with 
experiments^^; 6) universality with the correct scaling 
given by Tk- 

At the same time, due to the second order expansion 
of Eq. p3|). the theory does not take into account all 
inelastic cotunneling processes which are of the fourth 
order in the tunneling matrix elements. Therefore, the 
present theory underestimates the differential conduc- 
tance at voltages such that e\V\ > gfiB\H\ (see Fig. [3]). 
An additional consequence of this approximation is that 
at high magnetic fields the differential conductance as a 
function of the bias voltage is oversensitive with respect 
to the temperature (see Fig. 2]). These drawbacks can be 
eliminated by taking into account higher order terms in 
Eq. ([15]). This is technically more complicated and will 
be addressed in our future research. 
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